### Open Data set in .txt
coliga=read.table("coligacoes.txt", header=T)

attach(coliga)

### Plot comparing NEP and Electoral Coalitions by state
plot(contingente~nep, xlab="NEP", ylab="Tamanho da Coligação Eleitoral")
points(contingente[estado=='Rio Grande do Sul']~nep[estado=='Rio Grande do Sul'], col="black", pch=20)
points(contingente[estado=='Minas Gerais']~nep[estado=='Minas Gerais'], col="red", pch=20)
abline(lm(contingente[estado=='Rio Grande do Sul']~nep[estado=='Rio Grande do Sul']), col="black")
abline(lm(contingente[estado=='Minas Gerais']~nep[estado=='Minas Gerais']), col="red")
legend("topright", c("Minas Gerais","Rio Grande \n do Sul"), pch=16, col=c("red","black"))

detach(coliga)

### Open Data Set Gabinetes and Save as .txt
gab=read.table("gab.txt", header=T)
names(gab)

attach(gab)
estado=as.integer(estado)

### Plot comparing Fragmentação Partidária and Number of Partners by state
plot(parceiros~fragment, xlab="Fragmentação Partidária", ylab="Número de Parceiros na Coalizão")
points(parceiros[estado==1]~fragment[estado==1], col="red", pch=20)
points(parceiros[estado==2]~fragment[estado==2], col="black", pch=20)
abline(lm(parceiros[estado==1]~fragment[estado==1]), col="red")
abline(lm(parceiros[estado==2]~fragment[estado==2]), col="black")
legend("topleft", c("Minas Gerais","Rio Grande \n do Sul"), pch=16, col=c("red","black"))

### Plot comparing NEP and Tamanho da coalizão
a <- lm(cadeiras~nep+I(nep^2)); summary(a)
gf= function(x) -485.02 + 137.702*x -8.65*x^2
x=seq(5.5,10,.5)

plot(cadeiras~nep, xlab="NEP", ylab="Tamanho da Coalizão", main="Correlação entre NEP e Tamanho da Coalizão \n em MG e RS (1999-2006)", xlim=c(5,10), ylim=c(0,80))
points(cadeiras[estado==1]~nep[estado==1], col="red", pch=20)
points(cadeiras[estado==2]~nep[estado==2], col="black", pch=20)
lines(x, gf(x), lty=2)

### Plot Tempo da Coalizão and Size da Coalizão
plot(meses~cadeiras, xlab="Tamanho da Coalizão", ylab="Tempo de Duração da Coalizão")
points(meses[estado==1]~cadeiras[estado==1], col="red", pch=20)
points(meses[estado==2]~cadeiras[estado==2], col="black", pch=20)
abline(lm(meses[estado==1]~cadeiras[estado==1]), col="red")
abline(lm(meses[estado==2]~cadeiras[estado==2]), col="black")
legend("topleft", c("Minas Gerais","Rio Grande \n do Sul"), pch=16, col=c("red","black"))